SINGLE-PERIOD  STOCHASTIC  INVENTORY 
PROBLEMS  WITH  QUADRATIC  COSTS 


Moon  Ho  Song 


"HOX  LIBRARY 
-  POSTGRADUATE  SCHOOQ 
FREY,  CALIFORNIA 


NAVAL  POSTGRADUATE  SCHOOL 

Monterey,  California 


THESIS 


SINGLE-PERIOD  STOCHASTIC  INVENTORY 
PROBLEMS  WITH  QUADRATIC  COSTS 

by 

Moon  Ho  Song 


Thesis  Advisor 


G.  F.  Lindsay 


March  1974 


Approved  Ion.  pubtic  h.eJL<mi>z;  dl&tnxbLLtion  wnZimLtzd. 


Single-Period  Stochastic  Inventory 
Problems  with  Quadratic  Costs 

by 


Moon  Ho  Song 
Major,  Republic  of  Korea  Army 
B.S.,  Republic  of  Korea  Military  Academy,  1961 


Submitted  in  partial  fulfillment  of  the 
requirements  for  the  degree  of 


MASTER  OF  SCIENCE  IN  OPERATIONS  RESEARCH 


from  the 

NAVAL  POSTGRADUATE  SCHOOL 
March  1974 


DUDLEY  KNOX  L.BRARV 

TEeEY«  CAUFORN/A  93940 


ABSTRACT 

Single-period  inventory  problems  such  as  the  news- 
paper boy  problem  having  quadratic  cost  functions  for 
either  or  both  shortages  and  overages  are  examined  to 
determine  the  optimal  order  level  under  various  prin- 
ciples of  choice  such  as  minimum  expected  cost,  aspira- 
tion level,  and  minimax  regret.   Procedures  for  finding 
the  optimum  order  levels  are  developed  for  both  con- 
tinuous and  discrete  demand  patterns. 
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I.    INTRODUCTION 


In  any  decision  of  interest,  there  are  two  or  more 
alternative  courses  of  action  among  which  the  decision 
maker  must  choose.   A  principle  of  choice  then  indicates 
which  alternative  is  actually  to  be  selected.  [2]   In 
inventory  problems  we  are  concerned  with  making  optimal 
decisions  with  respect  to  an  inventory  system  and  in 
particular,  with  making  optimal  decisions  that  minimize 
the  total  cost  of  an  inventory  system.  [3] 

The  purpose  of  this  thesis  is  to  identify  and  relate 
certain  principles  of  choice  to  single-period  inventory 
problems  with  quadratic  cost  functions  for  either  or  both 
shortages  and  overages,  in  such  a  manner  that  the  optimal 
order  level  could  be  determined  under  each  principle  of 
choice.   Such  principles  of  choice  as  expectation,  aspira- 
tion level,  and  minimax  regret  will  be  examined  in  this 
paper . 

The  thesis  proceeds  by  first  reviewing  the  classic 
single  period  inventory  problem.   Procedures  for  finding 
the  optimal  order  levels  under  each  of  the  three  prin- 
ciples of  choice  when  the  cost  functions  for  both  short- 
ages and  overages  are  quadratic  are  addressed  in  Chapter 

III  for  both  the  continuous  and  discrete  cases.   Chapter 

IV  includes  consideration  of  problems  .having  (1)  a 
linear  cost  function  for  overages  and  a  quadratic  for 


shortages  and  (2)  a  constant  cost  function  for  shortages 
and  a  quadratic  for  overages.   In  Chapter  V  there  are 
conclusions  and  recommendations  for  further  work.   The 
philosophical  implications  of  inventory  theory  are  beyond 
the  author's  intentions  and  are  not  considered  here. 


II.   THE  CLASSIC  SINGLE-PERIOD  INVENTORY  PROBLEM 

This  chapter  is  concerned  with  single-period  inven- 
tory problems  where  the  demand  for  a  period  is  a  random 
variable  having  either  a  known  or  an  unknown  probability 
distribution,  and  costs  are  linear.   This  problem  has 
been  studied  extensively,  [1]   [2]  and  forms  the  basis 
for  the  work  presented  in  later  chapters. 

In  an  inventory  situation  in  which  items  are  ordered 

at  the  beginning  of  the  period,  let  C   be  the  unit  cost 

of  shortage  (or  profit  per  item  as  an  opportunity  cost) 

and  let  C   be  the  unit  cost  of  surplus  at  the  end  of  the 
s 

period.   The  decision  variable  is  S,  the  quantity  on 
hand  at  the  beginning  of  the  period.   Let  demand  D  be  a 
random  variable  which  denotes  the  demand  during  the 
period,  with  probability  density  function  f (D)  and 
distribution  function  F(D).   The  cost  per  period  is 


C(S)  = 


{C  (S-D) ,  0  <  D  <  S, 
-     - 
C  (D-S) ,  S  <  D, 


or  alternatively,  the  profit  per  period  is 


7T(S)  = 


C  D-C  (S-D) ,  0  <  D  <  S, 
s  —    — 


{C  D- 
CoS' 


S  <  D. 


Well  known  results  regarding  the  optimal  order  level 
S   may  be  summarized  as  follows. 


A.   OPTIMAL  SOLUTIONS  FOR  CASES  WHERE  F(D)  IS  ESTIMABLE 

To  minimize  expected  cost  or  maximize  expected  profit, 

the  optimal  order  level  S   is  chosen  so  that 

o 

F(S  -1)  <C/(C   +C)<F(S), 
o o    s     o        o 

if  the  demand  is  discrete,  or 

F(S  )  =  C  /(C   +  C  )  , 
o      o    s     o 

if  the  demand  is  continuous.  [1] 

To  maximize  the  probability  that  cost  does  not  exceed 

the  decision  maker's  aspiration  level  A,  where  the  demand 

is  continuous,  the  optimum  order  level  S   is  chosen  so 

o 

that  [2] 


f(S   -  A/C  )  =  f(S   +  A/C  )   . 
o       s        o       o 


B.    OPTIMAL  SOLUTIONS  FOR  CASES  WHERE  D      IS  ESTIMABLE 
BUT  F (D)  IS  NOT  maX 

To  minimax  regret,  the  optimal  order  level  S   is 
3  r  o 

chosen  such  that 


S   <C(D     +1)/(C   +C)<S   +1, 
o     o   max         so      o 


if  the  demand  is  discrete,  or 

S^  =  C  (D    )/(C   +  C  )   , 
o     o   max     s     o 

if  the  demand  is  continuous.  [2] 

The  next  two  chapters  will  be  deverloped  in  a  similar 
way  as  we  seek  solutions  for  the  quadratic  case. 


III.    SINGLE-PERIOD  INVENTORY  MODELS 
HAVING  QUADRATIC  COST  FUNCTIONS 

In  order  to  translate  a  realistic  inventory  problem 
into  a  mathematical  problem  of  minimizing  a  cost  function, 
both  flexible  and  simple  approximations  to  a  wide  range 
of  cost  relationships  are  desirable  to  allow  easy  mathe- 
matical solutions.   Consideration  of  the  kinds  of  costs 
involved  suggests  that  a  U-shaped  cost  curve  is  required. 
[6]   For  example,  the  cost  of  inventory  is  high  when  in- 
ventory is  large,  and  high  also  at  the  other  extreme  when 
inventory  is  so  small  that  there  are  frequent  runouts  of 
inventory.   Somewhere  between  these  extremes,  the  combined 
costs  are  at  a  minimum.   With  these  considerations  in  view, 
the  cost  functions  may  sometimes  be  approximated  with 
reasonable  accuracy  by  a  positive  definite  quadratic  form. 
[6] 

This  chapter  is  concerned  with  single-period  inventory 
models  having  quadratic  cost  functions  for  both  shortages 
and  overages  in  connection  with  both  continuous  and  dis- 
crete demand  patterns. 

A.    CONTINUOUS  DEMAND  PATTERNS 

The  procedure  developed  below  is  a  method  for  finding 

the  optimal  order  level  S   when  the  cost  functions  can  be 

o 

approximated  by 

-  D)    ,   0  <_    D  <_  S  , 
cost  per  period,  C(S)  =  {  (1) 

S)    ,   S  <  D  . 


(Cs(S  - 


1 .   Minimum  Expected  Cost  Solution 

Suppose  that  the  demand  for  a  period  is  a  random 

variable  having  a  known  probability  distribution.   It  is 

assumed  that  the  demand  is  continuous.   The  expected 

total  cost  per  period  of  the  system  is 

s  °° 

E(C(S))  =   /   C  (S-D)2  f(D)dD  +   /cq(D-S)2  f(D)dD.   (2) 

o  s 

To  find  the  optimal  order  level  S  ,  the  expected  cost 

r-  0 

function  is  .differentiated  with  respect  to  S  and  the 
results  are  set  equal  to  zero.   This  involves  differentia- 


B(t) 


tion  of  an  integral.   It  can  be  shown  that  if 

I 
F(t)  =   J  G(t,x) , 

f?  A(t) 

*  *? 
then 

B(t) 

8G(t,x)dx/8t+G[t,B(t)]dB(t)/dt  -  G  [t ,A (t) ] dA (t) /dt. 

A(t) 
This  is  known  as  Leibniz  Rule.   Applying  this  result  to 

equation  (2)  yields 

s  °° 

dE(C(S))/dS=2C   /   (S-D)f (D)dD  -  2C   /   (D-S)f (D)dD.  (3) 

o  s 

Setting  the  derivative  equal  to  zero  and  solving  leads  to 


(D  -S  )/D(S  )-S   F(S  )  =   C  -  C  /c  ,  (4) 

o      o    o     o        o    s'  o 


where  D  =  Expection  of  demand,  E(D),  and 


so 


D(Sq)  =  J       Df(D)dD.  (5) 

This  result  (4)  agrees  with  the  work  of  Sirichoke.  [5] 


To  check  whether  S   satisfying  equation  (4)  gives 
a  minimum,  the  second  derivative  must  be  examined.   From 
equation  (3 ) , 

s  oo 

d2E(C(S))/dS2  =  2Cs  /   f  (D)dD  +  2CQJ     f  (D)dD. 

O  s 

Since  f (D)  is  a  probability  density,  the  value  of  the 

second  derivative  is  always  non-negative  for  non-negative 

unit  costs  C   and  C  .   Thus  the  value  of  S   satisfying 
so  o        J 

equation  (4)  does  indeed  furnish  an  optimal  solution. 

Equation  (4)  is  not  tractable  in  general  because 
of  the  difficulties  imposed  by  D(S  ) .   For  the  uniform 
demand  distribution,  however,  it  is  possible  to  obtain  S 
and  the  minimum  cost  C   in  explicit  form.   This  is 
illustrated  as  follows.   Let  the  demand  density  be 


f(D)  = 


—    —   max 


{1/D      ,  0  <  D  <  D 
max      —    — 
0  ,       otherwise 


Then  D  =  D    /2  and 
max 


F  (S 


)=  f  f  (D)dD/D      =  S  /D 

o  J  max     o   max 


and 


D(S  )=  /    DrHBjdD/D 

o   J  max 


=  S   /2D 

o     max 


Substituting  these  values  in  equation  (4)  gives  the 
optimal  order  quantity 


S   =  D    /(l  +  -v/c  /C 
o     max        V  s   o 


)  . 


(5) 
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With  the  uniform  distribution  and  S  ,  equation  (2)  gives 

C=(CS    +C(D     -S  )3)/3D 

so      o   max     o       max 

as  the  minimum  expected  cost. 

2 .   Aspiration  Level  Solutions 

It  is  possibly  true  that  some  form  of  aspiration 
level  principle  is  the  most  widely  used  of  all  principles 
in  management  decision  making  as  alternatives  become  in- 
creasingly expensive  to  discover.   An  aspiration  level  is 
simply  some  level  of  cost  which  the  decision  maker  desires 
not  to  exceed.   For  the  inventory  problem  we  are  con- 
sidering, an  aspiration  level  policy  might  be  expressed 
as  follows.   For  a  given  aspiration  level,  A,  select  the 

optimal  order  level  S   which  maximizes  the  probability 

o 


Demand , D 


Figure  1.   Quadratic  Cost  Function  v.s.  Demand,  Showing 
Aspiration  Level 
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From  Figure  1,  it  follows  that 


pr(cost<    A)  =pr  (S-  Ja/C       <_    D    <_    S+  -v/a/C     )  = 


F(S+^A/Cq)-    F(S-^A/Cs     ) 


It  is  assumed  that  the  probability  density  func- 
tion of  the  demand  is  unimodal  with  respect  to  the  maxi- 
mum of  that  function  defined  on  the  range  of  demand.   To 

find  the  optimal  level  S  ,  pr(cost  <  A)  is  differentiated 

o  — 

with  respect  to  S  and  the  results  are  set  equal  to  zero: 


dpr(cost<    A)/DS    =    f(S+*/A/C     )     -f  (S-  */A/C     )     =    0     . 


This    leads    to 


f(S    +  VA/C    )     =    f(S    -  Va/c    ) 

O  O  OS 


as  our  basis  for  choosing  an  optimal  value  of  S. 

It  can  be  easily  shown  that  if  f(D)  is  unimodal 
and  symmetric  about  D  (as,  for  example,  the  normal  distri' 
bution) ,  then 


S   =  D  +  -(\/A/C   -Va/C  )  . 
o         2       s        o 


This  is  shown  in  Figure  2. 

If  f (D)  is  unimodal  with  mode  at  0  (as,  for 
example,  the  exponential  distribution),  then 


S   =  Va/c 
o         s 
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as  the  aspiration  level  solution.   This  is  shown  in 
Figure  3.   We  note  that  in  the  latter  case  it  is  unneces 
sary  to  estimate  the  unit  outage  cost  C   to  achieve  an 
optimal  solution. 


f  (D) 


>  Demand , D 


S   -  v/A/c 
o    v     s 


D  S  +  /A/C 
o  V     o 


Figure  2.   Application  of  the  Aspiration  Level  Decision 
Rule  to  a  Unimodal,  Symmetric  Probability 
Distribution  for  Demand 


f  (D) 


VVc 


Demand , D 


Figure  3.   Application  of  the  Aspiration  Level  Decision 
Rule  to  a  Unimodel  Probability  Distribution 
for  Demand  Having  Mode  at  0. 
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3 .   Minimax  Regret  Solutions 

A  decision  for  which  the  analyst  elects  to  con- 
sider several  possible  futures,  whose  probabilities  can- 
not be  estimated  is  called  a  decision  under  uncertainty. 
[2]   When  we  are  unable  to  estimate  f (D) ,  the  inventory 
problem  may  be  treated  as  a  decision  under  uncertainty. 
A  principle  of  choice  for  decisions  under  uncertainty 
has  been  proposed  by  L.  J.  Savage,  who  suggests  that  a 
new  matrix  called  a  "regret  matrix"  be  computed  first. 
For  each  possible  future  combination  of  demand  and  order 
level,  the  difference  should  be  computed  between  the 
actual  cost  that  will  occur  and  the  best  cost  that  could 
be  occurred  for  the  future  under  consideration.   This 
difference  is  called  "regret." [4] 

This  principle  of  choice  may  be  appropriate  for 

cases  where  the  maximum  demand,  D    ,  is  known  but  f (D) 

max 

is  not.   By  observing  the  cost  function  and  bearing  in 

mind  that  demand  takes  on  any  value  between  0  and  D    , 

max 

it  should  be  easily  seen  that  the  minimum  cost  will  be 
zero  for  any  given  order  level  S,  when  demand  turns  out 

to  be  the  order  level.   The  maximum  surplus  cost  will  be 

2 

C  S  ,  when  demand  is  zero.   The  maximum  shortage  cost 
s 

2 
will  be  C  (D    -  S)  ,  when  demand  turns  out  to  be  D 

o   max  max 

For  any  order  level  S  the  maximum  cost  will  simply  be 
the  maximum  of  these  two  quantities  and  so  will  the  maxi- 
mum regret.   It  will  be  noted  that  the  maximum  surplus 
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2 
cost  C  S   will  increase  as  order  level  increases  and  the 

s 

maximum  shortage  cost  will  decrease  as  order  level  in- 
creases.  The  minimax  cost  will  occur  when  they  are  equal, 
thus  S   satisfying  the  equation 

2  2 

C  S    =  C  (D     -  S  ) 
so      o   max     o 


will  be  the  optimal  order  level.   Solving  this  for  S   we 

o 

have 


S   =  D    /(l  +  v/C  /c  ) 
o     max       v  s   o 

as  the  minimax  regret  solution.   It  should  be  noted  that 
this  result  agrees  with  equation  (5)  which  gives  the  mini- 
mum expected  cost  solution  when  demand  is  uniformly  distri' 

buted  between  0  and  D    .   Accordingly,  for  the  case  where 

max 

f  (D)  cannot  be  estimated,  the  same  result  would  occur  if 
we  had  used  Laplace's  Principle  of  Insufficient  Reason 
which  suggests  computing  a  simple  average  when  probabili- 
ties cannot  be  estimated.  [2] 

B.    DISCRETE  DEMAND  PATTERNS 

Suppose  that  the  demand  is  a  discrete  random  variable 
having  a  known  probability  distribution.   When  demand  D 
and  order  level  S  are  constrained  to  discrete  units  0,  u, 
2u,  ...  and  p(D)  is  the  probability  mass  function  of 
demand,  then  the  expected  total  cost  of  the  system  becomes 

E(C(S))  =  C   X)  (S-D)  p(D)+  C   V    (D-S)  p(D).      (6) 
s  o  *— ' 

D=0  D=S+u 
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The  necessary  conditions  for  S   to  be  the  optimal  order 
level  are 


E (C(S  +  u) )  -  E (C(S  ) )  >  0 

o  o  — 

E(C(S  -  u) )  -  E (C (S  ) )  >  0  . 

o  o  — 


(7) 
(8) 


To  find  the  conditions  for  the  system,  tfaie  differences 


and 


E(g(S+u) )  -  E(C(S) ) 


E(C(S-u) )  -  E(C(S) ) 


must  be  evaluated  in  general  first.   Here, 

S 
E(C(S+u) )-E(C(S) )=2uC   £(S-D+  -u)p(D) 


D=0 


°°  1 

-2uC   £     (D-S-  yu)p(DD 


D=S  +  u 


and 


V^        1 

E(C(S-u)  )-E (C(S)  )=-2uC  2^  (S-D-  -u)p(D) 

SD=0        Z 

oo  , 

+  2uC   J2  (D-S+  -u)p(D). 

o  *"**  2 

D  =  S  +  u 

Applying  these  results  in  equations  (7)  and  (8)  and 
solving  leads  to 


L(S  )  <  (C   -  C  )/C     U(Srt)  , 
o   —    o     s    o  —     o 


(9) 


where 


L(S  )  =   (S  -  \u    -D)/(F(S  ) (S  -  iu)-D(S  )), 
o        o   2  o    o   2        o 

U(Sq)  =  (Sq+  ju  -5)/(F(SQ)  (Sq+  |-u)-D(So)), 
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and 


5(S  )  -  £   Dp(D)  - 


D=0 


It  can  be  easily  shown  that  these  are  also  sufficient 

conditions  for  S   to  be  the  optimal  order  level. 

o  c 

As  an  example,  suppose  u  =  5,  C   =  $3„  C   =  $30,  and 

5  O 

the  demand  for  the  item  for  a  period  is  random  with  proba- 
bilities p(0)  =  0.05,  p(5)  =  0.25,  p(10)  =  0.35,  p(15)  = 
0.30,  and  p(20)=  0  „  0  5  .   Computation  of  the  values  of  L(S) 
and  U(S)  are  summarized  in  Table  I.   Since 


(Co  *  Cs)/Co  =  °*9' 


the  optimal  order  quantity  is 


S   =  15, 
o 


because 


0.857  =  L(15)  <  0.9  =  (C  -C  ) /C   <  U(15)  =  0.98 

OS     o 


As  a  check  on  this  result,  E(C(S))  of  equation  (6)  is 
calculated  for  all  values  of  S.   This  is  done  in  the  last 
column  of  Table  I.   As  expected,  the  least  cost  of  the 
system  occurs  when  S  =  15. 

No  rule  is  given  for  an  aspiration  level  solution  in 
the  discrete  demand  case  because  it  is  straightforward 
to  compute  pr  (cost  <_  A)  for  each  order  quantity  S,  and 
choose  the  one  which  maximizes  pr  (cost  <  A) .   To 
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this,  we  will  consider  the  previous  example,   Let  A,  the 

aspiration  level  for  cost,  be  $1,000.   Then  the  optimal 

order  level  S   is  easily  found  as  15  units  from  Table  II 
o 

TABLE  I 
Tabulation  of  L(S)  and  U(S) 

D,S      p(D)      F(S)      D(S)      L(S)      U(S)      E  C(S) 


0 

0.05 

0.05 

0 

102 

-62 

$3 

,862.50 

5 

0.25 

0.30 

1.25 

15.5 

-2.75 

1 

,503.75 

10 

0.35 

0.65 

4.75 

-22 

0.668 

408.75 

15 

0.30 
0.05 

0.95 
1.00 

9.25 
10.  25 

0.857 
1.00 

0.98 
1.00 

172.  50 

20 

356.25 

TABLE  II 
Computation  of  prfcost  <  $1,000) 


D  0  5  10  15  20  pr(cost  <_ 

S   p(D)  0.05  0.25  0.35  0.30  0.05  $1,000) 

0  0  750  3000  6750  12000  0.3 

5  75  0  750  3000  6750  0.65 

10  300  75  0  750  3000  0.95 

15  675  300  75  0  750  1.00 

20  1200  675  300  75  0  0.95 


In  a  discrete  demand  pattern,  the  Savage  principle, 

which  suggests  computing  a  regret  matrix  first,  may  also 

be  appropriate  for  the  cases  where  D     is  estimable  but 

max 

p(D)  is  not.   For  each  possible  future  combination  of 
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demand  D  and  order  quantity  S,  the  difference  is  computed 
between  the  actual  cost  that  will  occur  and  the  minimum 
cost  that  could  be  occurred.   Having  completed  the  regret 
matrix,  the  optimal  order  level  is  selected  which  mini- 
mizes the  maximum  regret. 

To  illustrate  this,  we  will  again  consider  the  pre- 
vious example.   A  regret  matrix  is  completed  first  as 

shown  in  Table  III,  and  the  optimal  order  level  S   is 

c  o 

selected  as  15  units,  since  the  regret  of  that  quantity 
is  the  minimum. 

TABLE  III 
Regret  Matrix 


S  D 

0 

5 

10 

15 

20 

max  regret 

0 

0 

750 

3000 

6750 

12O00 

12000 

5 

75 

0 

750 

3000 

6750 

6750 

10 

300 

75 

0 

750 

3000 

3000 

15 

675 

300 

75 

0 

750 

750 

20 

1200 

675 

300 

75 

0 

1200 

So  far,  we  have  examined  the  cases  where  both  shortage 
and  overage  cost  functions  are  quadratic.   The  next  chap- 
ter will  be  devoted  to  examine  the  cases  where  cost  func- 
tions- are  either  a  combination  of  linear  and  squared  terms 
or  of  constant  and  squared  terms. 
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IV.    SINGLE  PERIOD  INVENTORY  MODEL.S  HAVING 
QUADRATIC  COST  FUNCTION  FOR  EITHER 
SHORTAGES  OR  OVERAGES 


Although  a  mathematical  form  for  describing  a  realis- 
tic inventory  problem  should  be  simple  and  easy  to  handle, 
it  must  be  flexible  enough  to  give  adequate  approximations 
to  a  wide  variety  of  situations.   With  these  considera- 
tions in  view,  the  cost  function  may  sometimes  be  approxi- 
mated by  either  a  combination  of  linear  and  squared  terms 
or  by  a  combination  of  constant  and  squared  terms  in  con- 
trolled and  uncontrolled  variables.  [6]    This  chapter  pro- 
ceeds by  first  considering  the  case  where  the  overage  cost 
function  is  linear  over  the  region  of  the  number  of  units 
overstocked  and  the  shortage  cost  may  be  approximated  by 
a  quadratic  function  of  the  number  of  units  short.   The 
case  where  the  shortage  cost  is  constant  and  the  average 
cost  is  quadratic  will  also  be  examined. 

A.    LINEAR  COSTS  FOR  OVERAGES  AND  QUADRATIC  FOR  SHORTAGES 
The  procedure  developed  below  is  a  method  for  finding 
the  optimal  order  level  S   when  the  cost  per  period  is 


C(S)  =  (C     (S-D)   ,  0  <_    D  <_  S, 


'  "  (10) 

^C  (D-S)   ,       S  <  D. 


1 .   Minimum  Expected  Cost  Solution 

When  the  demand  for  a  period  is  a  continuous 
random  variable  having  a  known  probability  distribution, 
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the  expected  total  cost  of  the  system  is 


"/ 


/ 


E(C(S))=    Cg (S-D) f (D)dD+  /  Cq(D-S)  f(D)dD. 
O  s 

Applying  Leibniz  Rule  to  equation  (11)  yields 


(ID 


dE (C (S) )/dS  =  C  F (S)  -  2C 

s 


./«- 


S) f (D)dD. 


(12) 


To  find  S  ,  equation  (12)  is  set  equal  to  zero 
leads  to 


This 


where 


C  /2C   =  (DC(S  )  -  S  FC(S  ) )/F(S  )   ,  (13) 

SO  o       o     o        o 


D°(So)  =  I       Df(D)dD 


and 


FC(S  ) 
o 


f 


f(D)  dD. 


To  check  whether  S   satisfying  equation  (13)  gives  a  mini- 
mum, the  second  derivative  must  be  examined.   From  equa- 
tion (12)  , 

oo 

d2E(C  (S)  )/dS2  =  C  f(S)  +2C   /  f(D)dD. 

s  o  J 

s 
Since  f (D)  is  a  probability  density,  the  value  of  the 

second  derivative  is  always  positive,  thus  the  order 

quantity  S   satisfying  equation  (13)  does  indeed  minimize 

E(C(S)  )  . 

Equation  (13)  is  not  tractable  in  many  cases  because 

—  Q 

of  the  difficulties  imposed  by  D  (S  ).   For  the  uniform 
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demand  distribution  it  is  possible  to  obtain  S   and  the 

o 

minimum  cost  C*  in  explicit  form.   This  is  illustrated 
as  follows.   Let  the  demand  density  be 


f(D)  = 


{1/D      ,  0  <  D  <  D      , 
max      —    —   max 
JO       ,  otherwise. 


Then 


and 


DC(S  )  =  -k)     -  (S  2/2D    )  r 
o  2    max      o     max 

F(S  )  =  S  /D 

o      o   max 


F**(S  )  «  1  -  (S  /D    )  . 

o  o   max 


Substituting  these  values  in  equation  (13)  gives  the 
optimal  order  quantity 


S       =    D  +(C    /2C     )-  \I{C    D  /C     )  +  (C    /2C     )2     .      (14) 

o  max  s         o         »       s    max       o  so 


The  minimum  cost  is 


C  =(C  S  2/2D    )+(C  (D    -S  ) ( (D    -S  )   (D    -S  )(D    -4S  ) 
s  o     max     o   max   o     max   o     max   o    max    o 


+  3S   ) )/3D 

o       max 


When  the  demand  for  a  period  is  a  discrete  random 
variable  having  a  known  probability  distribution,  the  ex- 
pected total  cost  of  the  system  is 


E(C(S)  )=C  Y,   (S-D)p(D)+C   S     (D-S)  p(D)   . 


(15) 


D=0 


D  =  S  +  1 
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Using  the  necessary  and  sufficient  conditions  (7)  and  (8) 
in  Chapter  III,  we  obtain  the  following  results: 


L(S  )  <  C  /2C   <  U(S  )  , 
o   —   s    o  —     o 


(16) 


where 


L(S  )  =  (DC(S  )+FC(S  ) (S  -  i) )/F(S  ) , 
o  o       o    o   2       o 

U(Sq)  =  (DC(So)+FC(So) (Sq+  I))/F(So) , 

OO 

DC(SQ)=  £         Dp(D) , 
D=S0+1 


and 


FC(So)=  £        p(D) . 


D=S0+1 
2 .   Aspiration  Level  Solutions 

Suppose  one  wishes  to  keep  the  cost  per  period 
less  than  or  equal  to  an  aspiration  level  A.   One  might 
want  to  choose  order  quantity  S  so  that  pr  (cost  <  A)  is 
maximized . 


Cost  per 
period ,  C  (S) 


S  -  A/C 


S  + 


yp^o 


Demand , D 


Figure  4.   Cost  v.s.  Demand  with  Linear  Costs  for  Overages 
and  Quadratic  Cost  for  Shortages,  Showing 
Aspiration  Level 
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From  Figure  4,  it  is  seen  that  pr (cost  <A)  is  equal  to 


»r  (  (S-(A/Cs)  )  <  D  <_  (S  +  <^A/Co)  ) 


and  thus 


pr(cost  <  A)=  F(S+-/a/C  )-  F(S-(A/C  )). 
—  \    o  s 


It  is  assumed  that  the  probability  density  function  of 
the  demand  is  unimodal  with  respect  to  the  maximum  of 
that  function  defined  on  the  range  of  demand.   To  find 
the  optimal  order  level  S  ,  pr  (cost  <  A)  is  differen- 
tiated with  respect  to  S.   Setting  the  result  equal  to 
zero  yields  the  rule 

f(S  +JA/C  )  =  f(S   -  (A/C  )). 
o  M  o        o        s 

From  this  it  can  be  easily  shown  that  if  f (D)  is  unimodal 
and  symmetric  about  D,  then 


S   =  D  +  -((A/C  )  -(/A/C  ) ) . 
O  2        S     ^      o 


Also,  if  f (D)  is  unimodal  with  mode  at  0,  then 


S   =  A/C   , 
o       s 


as  the  aspiration  level  solution. 

3 .   Minimax  Regret  Solutions 

Suppose  that  the  maximum  demand.  D    ,  is  esti- 
rr  max 

mable  but  f (D)  is  not  and  we  wish  to  minimize  the 


This  analysis  is  similar  to  that  of  Chapter  III. 
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maximum  regret.   By  the  same  argument  developed  in  Chap- 
ter III,  we  have 

max{c  S,  C  (D     -  S)2} 
s     o   max 

as  the  maximum  regret.   Since  the  maximum  surplus  cost 
C  S  increases  linearly  as  the  order  level  increases  and 
the  maximum  shortage  cost  decreases  as  the  order  level 
increases,  S   satisfying  the  equation 

C  S  =  C  (D     -  S  )   , 
s      o   max     o 


.» 


gives  the  optimal  order  level.   Solving  this  for  S   we 
^  o 

obtain 


S   =  D      +(C  /2C  )-  /(C  D /C  )+ (C  /2C  ) ^ 
o     max     so       s  max   o     so 


as  the  minimax  regret  solution.   This  result  agrees  with 
equation  ( 14 )  . 

If  the  demand  is  discrete  we  obtain  the  follow- 
ing result : 


S    -1    <D  +(C    /2C     )-     /(C     (D  +1)/C     )+(C    /2C     )2     <     S 

o         —   max         so  s       max  o  s         o         —      o 


as  the  optimal  order  level  solution. 


B.    CONSTANT  COST  FUNCTION  FOR  SHORTAGES  AND 
QUADRATIC  FOR  OVERAGES 

The  procedure  developed  below  is  a  method  for 

finding  the  optimal  order  S   when  the  cost  per  period 

is  (  C  (S-D)  2 ,     0  <_  D  <_  S, 

C(S)  =<   S  (17) 

K  ,  S  <  D, 
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where  the  cost  due  to  understock  at  the  end  of  the  period 
is  a  positive  constant  K,  regardless  of  the  number  short. 
Such  a  function  might  be  appropriate  for  cases  where 
shortages  cause  system  failure. 

1 .   Minimum  Expected  Cost  Solutions 

When  demand  for  a  period  is  a  continuous  random 
variable  having  a  known  probability  distribution  which 
is  dif f erentiable  everywhere  on  the  domain, 


E(C(S))=/c  (S-D) 2f (D)dD+  /  Kf(D)dD. 


(18) 


Here , 


dE(C(S))/dS  =  2C   /(S-D)f (D)dD-Kf (S 


(19) 


Setting  equation  (18)  equal  to  zero  and  solving  this 
leads  to 


K/2C   =  (S  F(S  )  -  D(S  ) )/f (S  ) , 
s      o    o        o       o 


(20) 


where 


D(S  )  =  /  Df (D)dD, 


and 


'."/ 


F(S_)  =      f(D)dD. 
o 

Since 


d2E(C(S) )/dS2  =  2C 


/ 


is  always  positive,  S   satisfying  equation  (20)  does  give 
the  optimal  solution. 
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For  the  uniform  demand  distribution  it  is  possible 

* 

to  obtain  S   and  the  minimum  cost  C   in  explicit  form. 

o  c 

This  is  illustrated  as  below,  where  the  demand  density  is 

f(D)     = 


{1/D  ,     0    <    D    <    D  , 

max  —        —      nriax 

0  ,    otherwise. 


Then 


and 


D(S  )=  S   /2D    , 
o     o     max 

F  (S  )=  S  /2D    , 
o     o    max 


f(S  )=  1/D 

o       max 


Substituting  these  values  in  equation  (20)  yields  the 
optimal  order  quantity 


S   =  /K/C   . 
o        s 


With  the  uniform  distribution  and  S  ,  equation  (18)  gives 

*  3 

C   =  (C  S    +  3K(D     -  S  ) )/3D 

s  o         max     o      max 

When  the  demand  for  a  period  is  a  discrete  random 
variable  having  a  known  probability  distribution  p(D), 
then 

S  CO 

E(C(S) )= £  Cg (S-D)2p(D)+  £     Kp(D). 
D=0  D=S+1 

Using  the  necessary  and  sufficient  conditions  (7)  and  (8) 

in  Chapter  III,  we  get 

L(S  )  <  K/2C   <  U(S  )  , 
o   —      s  —     o 
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where 


and 


L(S  )  =  (F(So)(Sq-  -)     -    D(So)l/p(SQ) f 


U(S  )  =  (F(S  )  (S  +  -)     -    D(S  ))/p(S  ) 
o         002        o       o 


2 .   Aspiration  Level  Solutions 

Our  cost  function  is  shown  in  Figure  5.   For  an 
aspiration  level  solution  we  will  be  concerned  with  the 
relative  magnitudes  of  constant  shortage  cost  K  and 
aspiration  level  A. 


S-/A  /C   S-/A  /C 


Demand , D 


Figure  5.   Cost  v.s.  Demand  with  Constant  Costs  for 

Shortages  and  Quadratic  for  Surplus,  Showing 
Aspiration  Level 
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When  A  >  K,  it  is  easily  seen  from  Figure  5  that 


pr(cost    <    A,)     =    pr(S-    /A,/C       <    D    <    °°) 
—      1  Is    — 


and  thus 


pr(cost  <_   A  )  =  1  -  F(S-  /A  /c  )  , 


and  the  probability  pr(cost  <_   A  )  is  maximized  when 


F(S-  /a. /C )  is  at  the  minimum,  which  will  occur  when 
1   s 

S  takes  any  value  such  that 


0  <  S   =  /A,/C   . 
—   o      Is 


When  A  <  K,  we  see  from  Figure  5  that 


pr(cost  <_  A  )  =  pr(S-  /A  /C   £  D  <_    S), 


or 


pr(cost  <  A  )  =  F(S)  -  F(S-  /A  /C  ) 


It  follows  that 


dpr(cost  <_   A2)/dS  =  f(S)-  f(S-  /a  /C  )  . 


Setting  the  result  equal  to  zero  yields  the  rule 


f (S  )  =  f (S   -  /A_/C  ) . 
o         o      2   s 


Thus  it  can  be  shown  that  if  f  (D)  is  unimodal  and  syni' 
metric  about  D  (e.g.,  normal),  then 


S^  =  D  +  -/A  /C   . 
o         2   2s 
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Also,  if  f (D)  is  unimodal  with  mode  at  0,  then 


S   =  /A„/C   , 
o      2'  s 

as  the  aspiration  level  solution.   This  result  agrees 
with  the  aspiration  level  solution  to  the  quadratic  sur- 
plus cost  function. 

3 .   Minimax  Regret  Solutions 

Suppose  that  D     is  estimable  but  f  (D)  is  not  and 
max 

we  wish  to  minimize  the  maximum  regret.   By  inspection  of 
the  cost  function  (17)  and  assuming  a  continuous  demand,  it 
is  easy  to  compute  the  maximum  costs  for  a  given  S  as 
shown  in  Table  IV,  and  it  will  be  noted  that  this  is 
equivalent  to  the  regret  matrix. 

TABLE  IV 
A  Computation  of  Costs  for  Various  Demands 

Demand, D    Surplus  Cost    Shortage  Cost    Max  Cost 

0  C  S2  0  C  S2 

s  s 

SO  0  0 

D>S  0  K  K 

From  Table  IV,  it  is  clear  that  the  maximum  regret  is 
simply 

max(c  S  ,  k)  . 
s 
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Since  the  maximum  surplus  cost  C  S   increases  with  S 

s 

while  the  maximum  shortage  cost  is  unchanged,  we  want 

S   such  that 
o 


C  S    <  K  , 

so   — 


and  solving  for  S   we  have 

o 


0  <  S   <  /K/C 
—   o  —      s 


as  the  optimal  order  level  solution. 

If  the  demand  is  discrete  we  will  minimax  regret 

2 
for  any  S  such  that  C  S   <  K,  thus  S   is  any  quantity 

satisfying  the  inequality 


0  <  S   <  /K/C   . 
—   o  —      s 


We  have  seen  so  far  how  the  principles  of  choice 
may  be  applied  to  single-period  inventory  problems  for 
decision  making  under  risk  and  under  uncertainty.   Along 
these  lines  some  suggestions  for  additional  work  are  made 
in  the  next  chapter. 
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V.   CONCLUSIONS  AND  RECOMMENDATIONS 
FOR  FURTHER  WORK 

We  have  reviewed  the  classic  single-period  inventory 
problem  in  the  second  chapter  and  in  a  similar  way,  pro- 
cedures for  finding  the  optimal  order  level  under  various 
principles  of  choice  have  been  developed  when  the  cost 
functions  are  quadratic  for  both  shortages  and  surplus  in 
Chapter  III.   In  Chapter  IV,  other  interesting  single- 
period  inventory  problems  have  been  examined. 

When  demand  is  a  continuous  random  variable,  Leibniz 

Rule  becomes  a  very  useful  tool  in  differentiating  of  an 

integral  of  which  the  integrand  is  quadratic  and  it  is 

noted  that  non-negative  demand  distribution  combined  with 

non-negative  shortage  and  surplus  costs  assures  the 

second  derivative  of  the  expected  cost  function  always 

positive  and  thus  the  minimum  expected  cost  is  easily 

found  by  simply  setting  the  first  derivative  equal  to 

zero  and  solving  for  the  optimal  order  level  S  .   When 

o 

demand  is  a  discrete  random  variable,  the  finite  differ- 
ence  inequations  become  very  useful  tools  to  find  the 
minimum  expected  cost  solution.   In  a  discrete  demand 
case,  no  rule  is  given  for  an  aspiration  level  solution, 
because  it  is  straightforward  to  compute  the  probability 
that  cost  is  less  than  or  equal  to  a  given  aspiration 
level  A,  for  each  order  quantity  S,  and  choose  the  one 
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which  maximizes  pr  (cost  <_  A)  .   In  a  continuous  demand 
case,  the  assumption  of  the  demand  density  being  a  uni- 
modal  is  a  crucial  thing  to  achieve  an  aspiration  level 
solution.   Also,  we  have  noted  that  the  general  rules 
for  the  minimum  expected  cost  solution  are  not  tractable 

in  many  cases  because  of  the  difficulties  imposed  by  the 

—         —  c 
integrals  D(S  )  or  D  (S  ).   We  have  seen,  however,  that 
o  o 

for  the  uniform  demand  distribution,  it  is  possible  to 
find  the  optimal  solution,  which  agrees  with  that  of  the 
minimax  regret,  when  the  maximum  demand  is  estimable. 
The  use  of  more  than  one  principle  of  choice  has  resulted 
in  several  possible  optimal  solutions. 

In  a  variety  of  actual  contexts  the  cost  function 
might  be  better  approximated  by  a  piece-wise  linear  or 
exponential  function.   It  must  also  be  pointed  out  that, 
in  decisions  under  risk,  not  only  the  expectation,  but 
perhaps  also  the  variance  and  other  parameters,  of  the 
distribution  should  be  taken  into  account.   For  example, 
if  two  different  order  levels  have  the  same  expected  cost, 
then  the  one  with  the  smaller  variance  of  the  expected 
cost  will  be  chosen.   It  is  recommended  that  further  work 
along  these  lines  be  conducted. 

This  thesis  has  examined  a  number  of  interesting  and 
possibly  realistic  extensions  of  the  classic  single-period 
inventory  model,  both  in  terms  of  variations  in  the  cost 
function  and  in  the  applications  of  various  principles  of 
choice.   It  is  hoped  that  this  work  will  be  useful  to  those 
interested  in  this  kind  of  inventory  problem. 
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